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In this paper, the relation between skin friction and heat transfer along windward 
sides of blunt-nosed bodies in hypersonic flows is investigated. The self-similar 
boundary layer analysis is accepted to figure out the distribution of the ratio of 
skin friction to heat transfer coefficients along the wall. It is theoretically obtained 
that the ratio depends linearly on the local slope angle of the wall surface, and an 
explicit analogy expression is presented for circular cylinders, although the linear 
distribution is also found for other nose shapes and even in gas flows with chemical 
reactions. Furthermore, based on the theoretical modelling of the second order shear 
and heat transfer terms in Burnett equations, a modified analogy is derived in the 
near continuum regime by considering the rarefied gas effects. And a bridge function 
is also constructed to describe the nonlinear analogy in the transition flow regime. 
At last, the direct simulation Monte Carlo method is used to validate the theoretical 
results. The general analogy, beyond the classical Reynolds analogy, is applicable to 
both flat plates and blunt-nosed bodies, in either continuous or rarefied hypersonic 
flows. 
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= Skin friction coefficient 2t w /p 00 U% 0 
= Proportional factor of the analogy relation, Eq. ((6)) 
= Heat transfer coefficient 2 q w /PooU^ 

= Non-dimensional velocity along x axis u/u e 
= Non-dimensional total enthalpy H/H e 
= Total enthalpy per unit mass ( J/kg) 

= Thermal conductivity ( W/{m ■ I\ )) 

= Mach number in free stream 
= Index of the power-law shape leading edge 
= Prandtl number 
= Wall pressure ( N/m 2 ) 

= Normalized pressure distribution function, Eq. m 


9—0 


Pressure gradient parameter at the stagnation point — d 
Heat flux to the wall ( W/m 2 ) 

Section radius of the body of revolution, Eq. CD ( m) 

Curvature radius of the surface, Fig. © ( m) 

Radius of the inscribed circle of the base of the body, Fig. © ( m) 
Free stream Reynolds number PooUooRc/Poo 
Local Reynolds number PooUooR/Poo 
Temperature (A) 

Stagnation temperature in free stream ( K ) 

Velocity in free stream (m/s) 

Flow velocity along x and y (m/s) 

Rarefied flow criterion M/// /Rei 

Local coordinate in physical space, Fig. © (to) 

Gradient parameter 2d(lnrt e )/<i(hi£) 

Specific heat ratio 
Coefficients in Eq. © 3 ) 

Local coordinate in Lees — Dorodnitsyn transformation space 
Slope angle of the surface, Fig. © 

Surface slope angle at the base of the body, Fig. © 

Viscosity (kg/(m ■ s)) 

Density (kg/m 3 ) 

Normalized heat transfer distribution function, Eq. m 
Skin shear stress (N/M 2 ) 

Index of the viscosity-temperature power law 


Subscripts 
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e 

free 
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oo 


= Stagnation point condition 
= Boundary layer edge condition 
= Free molecular flow condition 
= Reference state 
= Wall condition 
= Free stream condition 


Superscripts 

l 1 ) = First order approximation based on Navier-Storks-Fourier equations 

= Second order correction based on Burnett equations 


2 





I. Introduction 


There have long been concerns on the connection between skin friction and wall heat-transfer 
rate (or simply, heat transfer), as macroscopically the two quantities are related respectively to the 
normal gradient of velocity and temperature, and microscopically they are momentum and energy 
transports arising from molecular moves and collisions. The most famous result of this subject is 
the Reynolds analogy in the flat plate boundary layer flow problem, where the skin friction and 
heat transfer were found proportional along the surface. But such simple relationship does not 
exist for curved surfaces. For example, in flows past circular cylinders or spheres the heat transfer 
reaches its maximum at the stagnation point and diminishes monotonously downstream while the 
skin friction is zero at the stagnation point and varies non-monotonically downstream. With regard 
to the relation between skin friction and heat transfer, there has not yet been a theory suitable for 
the curved surfaces in either continuum or rarefied gas flows. 

Separately speaking, the heat transfer in hypersonic flows has received much more attentions 


than the skin friction due to early engineering requirements. Lees [1] first developed theories of the 
heat transfer to blunt-nosed bodies in hypersonic flows based on laminar boundary layer equations. 
The method made rational predictions in high Reynolds number flows. Later Fay and Riddell 
[2] gave more detailed discussions on the stagnation point heat transfer and developed prediction 
formulas with improved accuracy. Considering the downstream region of the stagnation point, 


Kemp et al. [fj made an improvement to extend the theory to more general conditions. In practice, 
empirical formulas of the heat transfer distribution were also constructed for typical nose shapes 
such as circular cylinders and spheres 0 , [hi]. 

Besides above boundary layer analyses in continuum flows, theoretical studies of heat transfer in 
rarefied gas flows have been carried out by different approaches. Cheng Q| accepted the thin viscous 
shock layer equations and obtained analytical expressions of the stagnation point heat transfer from 
the boundary layer flow to the free molecular flow. Wang et al. [j], [§] presented a theoretical 
modelling of the non-Fourier aeroheating performance at the stagnation point based on the Burnett 
equations. A control parameter W r was derived as the criterion of the local rarefied gas effects, and 
formulas based on the parameter were found to correlate the heat transfer both at the stagnation 
point and in the downstream region. 

Unlike the heat transfer, the skin friction is often neglected in continuum flows over large blunt 
bodies, for the friction drag is much less than the pressure drag in those conditions. Most of the 
existing studies concerns only about turbulence flows rather than laminar flows. However for wedges 
or cones with small angles in rarefied gas flows, the skin friction contributes a significant part, as 
much as 50%, in the total drag Q. Unfortunately, there is still no reliable theory we could use to 
predict the skin friction over curved surfaces. It will be meaningful if we find a general Reynolds 
analogy by using which we could estimate the skin friction based on available heat transfer prediction 
formulas, or vice versa. In Lees and the followers’ research on the aeroheating performance of blunt 
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bodies, the momentum equations were solved coupled with the energy equation of boundary layer 
flows, which offers a breakthrough point to analyze the relation between the skin friction and the 
heat transfer. 

In the present work, the ratio of skin friction to heat transfer along curved surfaces is firstly 
discussed based on the self-similar solution of boundary layer equations. An expression with simple 
form is obtained for circular cylinders as a typical example. Subsequently, an extended analogy is 
deduced in the near continuum flow regime by considering the non-linear shear and heat transfer in 
the Burnett approximation, and it is found the rarefied gas effects on the analogy are characterized 
by the rarefied flow criterion introduced in our previous study. As a preliminary study, the molecular 
vibration and chemical reaction effects are out of consideration in the theoretical analysis. The direct 
simulation Monte Carlo (DSMC) method [loj] is also used to simulate present flows to validate the 
theoretical results. 


II. Ratio of Skin Friction to Heat Transfer in Continuous Flows 
A. Linear Analogy for Blunt-nosed Bodies. 

Fig. [T| is a sketch illustrating the hypersonic flow over a blunt-nosed cylindrical body or body 
of revolution. The local coordinate is set on the wall. 



Fig. 1: Blunt-nosed body in a hypersonic flow. 


In order to seek the self-similar solutions of the boundary layer equations governing the 
flow around blunt-nosed bodies, Lees et al. [1] introduced the Lees — Dorodnitsyn coordinate 
transformation: 


£ — / Pe^ePe^O 
Jo 


rj = 


VMJo ° 

and the normalizations of velocity and temperature: 


rtfpdy 


(1) 


f'(v) = u/u e 
g(v) = H/H e . 


( 2 ) 
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where m = 0 for planar bodies and m = 1 for bodies of revolution. With the transformation, the 
boundary layer equations were simplified and the self-similar solution have been obtained under 
certain conditions [di [n1|. 

Defining the coefficients of skin friction and heat transfer as C/ = 2t w /P ooU^ and 
Ch = 2q w /p 00 U% 0 , respectively, and with the approximation H e « 4 U in the hypersonic limit 
(Mqo 1), we have: 


9 l = 2Pr ^lL 

Ch CoqPw 


( 3 ) 


from Eqs. m and ©. 

From the compressible Bernoulli’s equation, the streamwise velocity u e along the edge of the 
boundary layer depends mainly on the wall pressure p 
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( 4 ) 


On the basis of the symmetry and smoothness of the pressure distribution near the stagnation point, 
Eq. ([4]) can be linearized to: 


Ue 

U ." 


Po^r-^- ' 9 


27 


( 5 ) 


with Mqo 1 and 9 -C 1, where po = 


d jp/po) 

W 2 


. In fact it was found that the linear variation 


J 61=0 


of u e with 6 is valid downstream the stagnation point till 9 = Air/9 from Korobkin’s experimental 


data [lj, 113. Therefore in following analyses, Eq. ([5]) is accepted not only near the stagnation 
point, but also in the downstream region. 


Cohen and Reshotko [14] numerically solved the self-similar boundary layer equations and 
calculated 2f”/g' w . The value depends on T w /Tq and varies with /3 = 2d(lnu e )/o?(ln^) along the 


surface. Actually according to the calculations from Kemp et. al [3], g' w can be taken as a constant 
along the isothermal wall with a slowly changing curvature radius. And from Cohen and Reshotko’s 
solutions [l^| , the slight variation of /" is of the same order as that of g ' w , and thus it is reasonable 
to assume a constant 2 f^/g' w as long as the self-similar assumption is satisfied. 

As a result, the expression of Cf /Ch in Eq. ([5]) can be written as: 

Cf 


with the coefficient 


C h 


C r = 2 Pr 


= C r 9 


( 6 ) 


f" 

J U. 
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Po¬ 


ll 


being independent of the location. The equation indicates that for a variety of nose shapes in 
hypersonic flows, Cf /Ch is proportional to 9 along the windward surface as long as T w is a constant. 
Thus, we have a more general form of the analogy relation between skin friction and heat transfer 
which is not restricted to the flat plate flow, compared with the classical Reynolds analogy. 
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B. Analogy for Circular Cylinders. 

In order to further calculate the coefficient C r , the shape dependent f^/g' w and po should be 
specified. As a typical example, the corresponding values for circular cylinders are to be presented 
in the following. 

First, the pressure distribution can be calculated by the Newtonian-Busemann theory [11], and 
at the stagnation point, we have 


Po = 3 


( 7 ) 


Second, since 2 f!^/g' w is regarded as a constant along the surface, we will calculate the value at 
the stagnation point where (3 = 1 for a planar body. Then from Cohen and Reshotko’s numerical 
calculations fl4|, 2/" /g' w was obtained within 0 < T w /Tq < 1.2 as plotted in Fig. ©. Obviously 
the points with T w /Tq < 1 fall into a straight line, and a linear correlation 


2(1 - T w /To)f"/ g' w = 2.48 + 1.88^/To (8) 

fits well with the points. With Eqs. (0 and ( 0 ) applied to Eq. ©, we have 

CllCh , Pr f ^^^ 9 



T w /To 


Fig. 2: Linear correlation of 2(1 — T w /T$) f!^ /g' w at the stagnation point of two-dimensional 
(/? = 1) bodies against T w /To based on data in [14], 

The explicit relation in Eq. 0 gives a convenient formula to predict C//Ch around a circular 
cylinder in hypersonic flows. Comparisons are given in the following with DSMC numerical 
simulations of the nitrogen gas flows past blunt-nosed bodies in hypersonic speeds under = 
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5 ~ 25, Re oo = 300 ~ 1500, T w /Tq = 0.01 ~ 0.4. The simulations are carried out by using the 
source code as described in |8|. The molecular vibration and chemical reaction effects are excluded 
to correspond with the derivations in this paper. 



Moo 5 8 10 15 20 25 


Fig. 3: Distributions of Cf/Ch for circular cylinders under different Mach numbers with 
Tw/Tq = 0.03. 



T w /T 0 0.01 0.05 0.1 0.134 0.2 0.4 

Fig. 4: Distributions of Cf/Ch for circular cylinders in different wall temperatures with = 10, 
except that T w /T$ = 0.134 case is from Santos [l]] under M ^ = 12 and Re x = 214. 

The numerically computed distributions of Cf/Ch along windward surfaces of circular cylinders 
are presented in Fig. m and dH). In Fig. (j3j) , the numerical results from different all fit 

well with Eq. ©. showing a Mach number independence, except that the ratio with = 5, 
the lower limit of hypersonic flows, is slightly lower. Meantime cases with T w /Tq = 0.01 ~ 0.4 in 
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Fig. Q, as well as a reference result from Santos |9|j with T w /Tq = 0.134, indicate that Eq. ® 
precisely describes the variation of Cf/Ch with T w /Tq. Other flow parameters in simulations are 
R c = 2 x 10 _ 3 m,= 300/\,poo = 1.5 x 10 _ 3 fcp/m 3 , and Pr = 0.71, 7 = 1.4 are accepted in Eq. 


© for comparisons. 

The Reynolds analogy for the flat plate flow shows that Cf/Ch distributes uniformly along the 
surface as mentioned. Throughout the zero-thickness leading edge, the ratio can be expressed as: 

( 10 ) 


Cl 

c h 


7+1 


7 +1 T,„ 


27(7 - 1 )Ml 0 4 7 T 0 


The equation was derived from the free molecular theory, and with M 00 1, T w /Tq -C 1, its range 

of application could be extended from the free molecular flow to the boundary layer flow |ljj. 

In practice, a finite thickness and a blunt nose always exist in the leading edge of a flat plate 
as illustrated in Fig. ®. Assuming a cylindrically blunted nose in the front of the plate, the ratio 
Cf/Ch at the flat segment can be obtained by taking 9 = 7 t/2 in Eq. © as: 


C f _ vrPr /3( 7 - 1) 2.48 + 1.88T W /T 0 

C h 2 Y 27 1 — T w /Tq 1 j 

Values of Cf/Ch in Eq. (fill) and for zero-thickness flat plate are plotted in Fig. (J3J) against T w /Tq. 
The difference between them is less than 5% with 0 < T w /Tq < 0.1. This result indicates that the 
present linear analogy in flows around curved bodies is in fact consistent with the classical Reynolds 
analogy for flat plate flows. 



Fig. 5: Comparison between the classical Reynolds analogy and the linear analogy for curved 
surfaces. M^ = 10, 7 = 1.4 and Pr = 0.71 in Eq. uni) and CD- 


C. Analogy for other nose shapes and flow conditions. 

Although lack of explicit expressions at present, due to the mathematical complexity, the 
existence of the self-similar boundary layer flow could also be found near the wall surfaces of many 
other nose shapes as long as the variation of the radius of curvature is slow Hj. Flows over two 


dimensional wedges with a variety of shapes under = 10, T w /Tq = 0.05 are also simulated 




















to extend our discussions, as shown in Fig. @. Among these results, the shape of the power-law 
leading edge is expressed as y = cx n , and 9 C = 80° in the present simulations. In simulations 
Re oo = 200, R c = 2 x 10 -4 m,T 00 = 300/i, = 4 x 10 -3 fcg/m 3 . Distributions of Cf/Ch along 

surfaces of power-law shapes with n = 0.5, 0.7 are both linear with 9 as expected. The lines diverge 
slightly from each other but still fit with Eq. ([9]) generally. The discrepancies may be caused mainly 
by the variation of the wall pressure distributions. 

The present theory is not suitable for the rectangle cylinder with a round shoulder, since the 
boundary layer around the body is highly non-selfsimilar js]. As presented in Fig. ©, Cf/Ch 
increases in the vertical windward segment with 9 maintaining zero. A similar variation can be 
observed near the stagnation point of the power-law shape with n = 0.3, where the surface is also 
nearly vertical and Cf/Ch increases sharply. However, the linear variation of Cf/Ch with 9 is still 
observed in the round shoulder segment of the rectangle cylinder and the downstream region of the 
n = 0.3 power-law shape. This phenomenon indicates a possibility that the present linear analogy 
relation could be extended to wider situations. 



Fig. 6: Cf/Ch from DSMC simulations of different shapes of blunt-nosed bodies with = 10 , 
Re oo = 200 and T w /Tq = 0.05. 


Further more, Cf/Ch is also found proportional to 8 when considering the real gas effects. 


Blottner 


15j calculated the boundary layer equations for the equilibrium air flow over a hyperboloid 


under M*, = 20. Holman and Boyd (16] computed the dissociating air flow over a sphere under 
Moo = 25 and Kn^ = 0.002 (Krioo = \ 00 /R 1 where Aoo is the mean free path of molecules in the 
free stream) with both DSMC method and Navier-Storks(N-S) equations. Results of Cf/Ch from 
both hyperboloid and sphere cases are demonstrated in Fig. 0. Although the slopes are much 
different from that of the above-mentioned calorically perfect gas condition, the linear feature can 
be still clearly observed. 

Despite assumptions and simplifications, the linear distribution of Cf/Ch on the windward of 
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Fig. 7: Cf/Ch from numerical simulations with real gas effects. 


curved surfaces is revealed by theoretical analyses and shown by DSMC simulations with different 
types of shapes. For circular cylinder cases, with the available explicit f!/,/g' w and po , good 
agreements are observed between the present analytical prediction and numerical results. 


III. Ratio of Skin Friction to Heat Transfer in Rarefied Flows 

The analyses in the above section are based on the boundary layer assumption, and thus the 
theory is valid only for continuous flows. When the flow deviates from the continuum regime, the 
nonequilibrium of molecular collisions causes non-linear shear and heat transfer, and then the linear 
Newtonian shear and Fourier heat transfer in the N-S equations fail. If the deviation is small, the 
constitutive relation_can be corrected by bringing in the second order shear and heat transfer in 


Burnett equations [17] as has been suggested by Wang et al. [7|, |8j] . In order to obtain a more 
general analogy relation covering the rarefied flows, the second order shear stress and heat transfer 
will be studied firstly in the near continuum regime, and then, based on numerical validation and 
calibration, the results will be extended to the more rarefied flow regime to calculate Cf/Ch■ In 
fact, an early exploration has been carried out in our previous work on the flat plates leading edge 
flow problem [18j . 


Instead of directly solving the Burnett Equations, the second order effects could be evaluated 
based on a perturbation point of view, i.e. to use the non-linear constitutive relations in Burnett 
equations to analyze the flow field features predicted by the first order approximation for instance 
the boundary layer theory or the computation method of N-S equations. The original form of the 
Burnett equations can be found in Chapman and Cowling’s derivations [17]. For planar bodies, 


taking the assumptions dp/dy « 0, v\ w = 0, -4§n\ = — fv 


dv I 


dy I 


7~1 TOp, 


7 TqR 


m and the constant 


wall temperature T\ w = T w . Besides, in the near continuum regime we have 


du I 
dy I tu 


u w /R and 
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dT | 
Tdy I 


1/P, then the Burnett shear and heat transfer near a planar curved surface become: 


r( 2 ) 


d 2 T 


w 3 


p 2 dp dT 


pRT dddy 
2 3 pRpdOdy 


1 

+ 2 U *'pRpT dd dy 


1 p 2 d 2 u 
~ 2 A ~pR~d6dy 


— w 2 


/i 2 du du 


+ (20 2 - 30 5 ) 


p Pc>0 dp 

p 2 dTdu 


pTP dy dd 


( 12 ) 


As has been testified in flows past the leading edge of flat plates, the above simplifications, although 
not strictly, still reflect the essential features of non-linear terms [18] . 

The flow held predicted by the first order continuum theory and the gradients near the wall 
need to be given before using Eq. m for further discussions. Assuming the normalized pressure 
and heat transfer distribution functions as 

p{9) =p/p 0 s(0) = ( 13 ) 

where po and q w $ are the pressure and the heat flux at the stagnation point, respectively. Both p 
and c are functions of 9 and are able to be obtained with different approaches such as analytical 


theories or data correlations [1, 
be transformed to 


. Taking po ~ tM^Poo in hypersonic limit, then Eq. in can 


dT 

'yM 00 pooP -gjj 


h (1) 

V10,0 


And with the linear analogy equation |6]), the normal gradient of u becomes 


du 

dy 


,kdT\ C ' 
pdy ) W U 00 


(14) 


(15) 


With Eqs. (fl4ll and (fl5l) submitted, Eq. (fl2ll becomes: 
( 2 ) _ + 0.5o;4pe/p / p 2 Uoo du 


9C r 


\kpRT dy 


) - 2w 2 VV(7-1) 


/ p 2 du 
\pUooR dy 


« § [0 3 C r (9p e /p) - 0 4 CV(1 + 0?®/?) + 4Pr(20 2 - 30 5 ) V(7 - 1)Ay] 




pRU a 


, dT 

: k si 


(16) 


where qg = dc/d0, pe = dp/d9. In Eq. (flfil) . expressions of the second order shear stress and heat 
transfer contain p^ and k^ explicitly, and thus, more meaningfully, we could get the relative 
magnitudes of the second order effects 


r (2) 

1 W 

7 11 

1 w 

gw' 1 

qw' 1 


^A + 0-5 U4P9/P p 2 r Uoo _ / n ^ Hr 

KprRTr PrUooR 


9C r 


(17) 


[0iC r (l + 9qe/q) ~ 9 3 C r (9p g /p) - 4Pr(20 2 - 30 5 )\/(7 - 1)/^] ■ 


\lr 


p r RU oo 

The subscript r represents the quantities under the reference temperature T r , as in hypersonic flows 
the reference temperature method is usually employed to present characteristics of boundary layers. 
Here T r = \{T W + To), and the viscosity-temperature power law p r /Poo = {T r /T 0 0 ) u is applied. 
Then Eq. (flTl) becomes 


Tw' 1 
Tw ^ 
qw^ 
qw ) 


= T 1 (9)W r 
= r 2 (9)w r 


(18) 
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where W r = M^/Rei, Rei = PooRUoo/ Moo, Ti and T 2 are: 


r i = 


r 2 = 


Pr {7 - l) 1+ “ + 0.5w 4 i5e/p / T w + T 0 


7P 


6»a 


4Tn 


— 2^2 


_ 1 W+l/2 


(7- 1) 


yl/2 


(7 2 ^- 1+ " [ 04^(1 + **) - e 3 C r e -f - APr(26 2 - 30 s )V(7-l)/7] ( 


T w + To 

4T 0 

Tw + To 

4T n 


(D+l 


CJ + 1 


(19) 


With ri, and being added to calculate the skin friction and heat transfer, we can introduce a 
modification to the linear analogy Eq. ([ 6 ]): 


Cl r^_ 

Ch 00 & 


J2) 


1 +T 1 W r 


.( 2 ) 

Qw 


C r 0 


( 20 ) 


1 + r 2 w r 

Values of Fi and T 2 denote the influence of the body shape and vary with 0, and generally both 
of them are in the order of unit. For a typical case of the hypersonic nitrogen gas flows past a 
circular cylinder under T w /Tq ks 0, with the modified Newtonian pressure and the heat transfer 
fitting formula from Beckwith and Gallagher [ 3 ] accepted in Eq. (d, an approximation can be 
taken in the range ofO<0<7r/6as 


Ei — r 2 » 0.43 + O.350 2 + O.180 4 


from the Taylor series expansion. Then in the near continuum regime with W r 1, Eq. m can 
be simplified to 

% « [1 + (Er - r 2 )W r ]C r e (21) 

The first order correction (Ti — E 2 )lE r is similar with Wang et al.’s result in the stagnation 
point heat transfer problem Q. It turns out that, W r is a control parameter of the rarefied gas 
effects not only at the stagnation point but also in the downstream region. Although the correction 
factor in Eq. (l2lT) is explicit and clear, it will lose its credibility when the rarefaction degree of 
the flow is sufficiently high, and due to the complexity of the transition flow, the data fitting and 
calibrating seem still unavoidable to get a practical general analogy. In fact, based on the rarefied 
flow criterion W r , a bridge function can be built between the continuum limit Eq. © and the free 
molecular limit (C / /Ch)free ~ 2 sinO. For circular cylinders, the function is carried out as 

Cf C r 6 + lAsindW r . 

C h l + 0.7W r y ’ 

Eq. EH with different W r is plotted in Fig. © compared with the DSMC simulations of nitrogen 
gas flows past circular cylinders. It can be seen that the variation of Cf/Ch with 6 is no longer 
linear in the rarefied gas flow regime. 


IV. Conclusions and Remarks 

The relation between skin friction and heat transfer for blunt-nosed bodies in hypersonic flows, 
named the general Reynolds analogy, has been investigated in this paper by using the theoretical 
modelling and the DSMC methods. 
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Wr 0.2 0.5 2 5 20 

Fig. 8: Non-linear analogy in hypersonic rarefied flows past circular cylinders with Moo = 10 ~ 25, 
Tw/Tq = 0 . 01 . 


First, based on the boundary layer flow properties, the ratio of the skin friction to the heat 
transfer for the blunt-nosed body was found proportional to the local surface slope angle. As a 
typical demonstration, an explicit expression of the ratio was derived for circular cylinders. Also, 
numerical calculations indicated that this characteristic exists for other blunt-nosed shapes even in 
chemically reactive flows. 

Second, the analogy in rarefied gas flows was analyzed. In the rarefied flow regime, the deviation 
from the linear distribution of the ratio was proved to be controlled by the rarefied flow criterion W r . 
Therefore, a bridge function was constructed based on W r to describe the analogy in the transition 
flow regime. 

This study, combined with our former investigation on the flat plate leading edge flows [jjj], 
clarifies the general Reynolds analogy in the whole flow regime for both flat plates and blunt-nosed 
bodies in hypersonic flows. 

The present general Reynolds analogy has potential usefulness in the engineering practice. From 
the analogy relation, the skin friction is related to the heat flux along surfaces, or further to the 
heat flux at the stagnation point and its normalized distribution downstream. As a result, the 
viscous drag integrated from the skin friction is proportional to the stagnation point heat flux, 
which suggests that if one of them is known, the other could also be obtained immediately. 
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